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ABSTRACT 

Thermonuclear burning on the surface of a neutron star causes the expansion of a thin outer layer 
of the star, AR(t). The layer rotates slower than the star due to angular momentum conservation. 
The shear between the star and the layer acts to twist the star's dipole magnetic field giving at first 
a trailing spiral field. The twist of the field acts in turn to 'torque up' the layer increasing its specific 
angular momentum. As the layer cools and contracts, its excess specific angular momentum causes it 
to rotate faster than the star which gives a leading spiral magnetic field. The process repeats, giving 
rise to torsional oscillations. We derive equations for the angular velocity and magnetic field of the layer 
taking into account the diffusivity and viscosity which are probably due to turbulence. The magnetic 
field causes a nonuniformity of the star's photosphere (at the top of the heated layer), and this gives rise 
to the observed X-ray oscillations. The fact that the layer periodically rotates faster than the star means 
that the X-ray oscillation frequency may "overshoot" the star's rotation frequency. Comparison of the 
theory is made with observations of Chakrabarty et al. (2003) of an X-ray burst of SAX J1804. 4-3658. 

Subject headings: neutron stars: magnetic fields: type I X-ray bursts 
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1. INTRODUCTION 

Type I X-ray bursts have been observed in a number 
of low-mass X-ray binary systems. They are character- 
ized by a rapid (~ 1 — 10 s) rise in the observed flux, 
followed by a relatively slow (~ 10 — 100 s) decline (see re- 
views by Lewin et al. 1995; Strohmayer & Bildsten 2003). 
Some objects show highly coherent oscillations during the 
burst. The oscillation frequency varies slightly during the 
burst, but asymptotically approaches the neutron star ro- 
tation frequency (Strohmayer & Bildsten 2003 and refer- 
ences therein; van der Klis 2000). The oscillations indicate 
nonuniform emission from the star's photosphere, surface 
of the star. The nonuniformity is likely due to the star's 
non-aligned dipole magnetic field. 

It is generally thought that Type I X-ray bursts are 
caused by thermonuclear explosions on the surface of a 
neutron star. If the neutron star's magnetic field is strong 
enough to channel accretion onto the star's surface, then 
the thermonuclear burning will not be uniform over the 
surface. It is then likely that the oscillations are produced 
by spin modulation of the burst flux. The exact mecha- 
nism of the oscillations, however, remains unclear in spite 
of a number of studies (Cumming et al. 2002; Spitkovsky 
et al. 2002; Bhattacharya et al. 2005; Piro & Bildsten 
2005). A significant difficulty encountered by models is 
explaining the oscillation frequency drifts. Cumming et al. 
(2002) investigated the effect of hydrostatic expansion of 
the neutron star atmosphere on the burst oscillations. An 
expanded layer will rotate slower than the star in order 
to conserve angular momentum, and then return to the 
stellar rotation rate during the contraction phase. This 
could cause the drift in the oscillation frequency. How- 
ever, Cumming et al. concluded that this phenomenon 
alone is not enough to explain the largest observed fre- 
quency drifts. They also noted that a strong stellar mag- 
netic field could inhibit angular velocity changes in the at- 



mosphere. In some bursts the oscillation frequency 'over- 
shoots' its asymptotic frequency (Strohmayer 1999; Miller 
2000; Chakrabarty et al. 2003). 

Here, we investigate the influence of a star's magnetic 
field on the rotation of its heated outer layer. The theory is 
developed in §2 where the induction equation is combined 
with angular momentum conservation to obtain equations 
for the rotation rate and magnetic field in the heated layer. 
Later in §2 we include the influence of the magnetic diffu- 
sivity and the viscosity of the heated layer. In §3 we discuss 
sample solutions for the rotation rate and magnetic field 
in the layer for a case relevant to a burst of SAX J1808.4- 
3658 described by Chakrabarty et al. (2003). Section 4 
gives the conclusions of this work. 

2. THEORY 

We consider that a thermonuclear explosion on the sur- 
face of the star creates a heated layer of mass AM = const 
and thickness AR(t) outside of the star's equilibrium ra- 
dius with AR considered for simplicity to be uniform 
over the star's surface. The actual layer may be an equato- 
rial band (Inogamov & Sunyaev 1999). The heating of the 
layer occurs rapidly and it expands rapidly giving AR(0). 
Subsequently, the layer slowly cools and AR(t) evolves and 
then decreases on a time-scale 10—100 s. We consider that 
AR(t) is a known function of time obtained from the ra- 
diative transfer and the radial force balance. 

We use a spherical inertial coordinate system (R, 0, 
and a coordinate system (R,0,cj>') rotating with the star 
at the constant angular rate £7* . The flow field in the layer 
is 

v = vrR + v^4> ■ (1) 
For the radial velocity we assume 
R — i?* 



AR (t) 
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(2) 



where AR= dAR/dt. The azimuthal velocity can be writ- 
ten as 



(3) 



where v'^R*, 9) = 0. 

We consider the case of an orthogonal rotator where the 
magnetic moment fi is perpendicular to the rotation axis 
fl* as sketched in Figure 1. This situation would give X- 
ray oscillations at about twice the rotation frequency of 
the star. This case is valuable in that it is amenable to 
an analytic treatment, and it indicates the behavior for 
non-orthogonal rotators where the X-ray oscillations are 
at about the rotation frequency of the star. For a non- 
orthogonal rotator, we suggest that the value of /i in the 
following expressions be replaced by fj,± = \fi x f2\/\S2\. 

We focus on the equatorial region of the star's surface, 
\6 — n/2\ <7r/6. The evolution of the magnetic field within 
the layer is described by the induction equation which fol- 
lows from Faraday's law and infinite conductivity, 




~dt 



= Vx(v x B) 



(4) 



with v given by equation (1). Later, we include a finite 
conductivity. Inside the star R < i?» the magnetic field 
is frozen-in and rotates with the star. Outside the layer, 
R > R« + AR, the magnetic field is considered to be a 
vacuum field which rotates but is otherwise unaffected by 
processes in the layer. Any field lines linking the star and 
disk are expected to be opened owing the large difference 
between the angular velocity of the inner disk and the star 
(Lovelace, Romanova, & Bisnovatyi-Kogan 1995). 

It is useful to describe the magnetic field in terms of a 
flux function ^ = RAg, where A = (0, Ag, 0) is the vector 
potential. Therefore, in the equatorial region of the star, 



B R = - -t^t 



i a* 

R^~d^ 



B, 



-— (5) 
RdR' (5) 



and Bg = 0. We adopt a reference frame rotating with 
the star with <f>' = <j> — Q*t the azimuth relative to a fixed 
point on the star. In place of equation (4) we have 



= -(v • V)* = -vn 

dt 1 ; OR Rd<j>' 



(6) 



where v' = (vr, 0,v'^) and v'^ = — 17* i?*. This equation 
simply says that the flux function is advected with the 
flow. 



Fig. 1. — The bottom panel shows the initial configuration and 
the top panel shows the configuration after some time. AR(t) is the 
thickness of the heated layer, and N and S indicate the north and 
south magnetic poles of the star. The vacuum field outside of the 
layer (R > R„ + AR) rotates but does not change its form. 



We solve equation (6) by taking 



y(R 7 <t>',t)=iB Rlexp 



i / dR'k R {R',t)+i<j>' 



R, 



(7) 



The physical solution is the real part of ^, denoted 3?(*); 
B is a constant field strength at the magnetic poles of the 
star; and kn is the radial wavenumber which remains to 
be determined. This form of solution was proposed earlier 
for magnetized disks (Bisnovatyi-Kogan & Lovelace 2000). 
From equation (5) we have 

Br = B § exp(...) , B = exp(...) , (8) 

where the exponential factors are the same as in equation 
(7). Initially, there is no toroidal magnetic field so that 

M* = o) = o. 

Substitution of equation (7) into (6) gives 



-k R 



R — -R* 



ir, dt A R 
Taking the R— derivative of this equation gives 



(9) 



dk 



R 



kj> AR R- R* 



AR 



dk 



R 



dt AR AR dR dR\R 

In the following we assume that \ {AR/kR)(dkR/ dR) \ <C 1. 
As a result the middle term of the right-hand side of this 
equation can be neglected. 

Note that Su(R, t) = v'^/R is the difference between the 
angular velocity of the layer and the star. We let t be 
the specific angular momentum of the layer matter and 
1* = Rl^l* the specific angular momentum of matter in 
the equatorial region of the star. We then have 

5i = e-e* = 2(r- R*)R*n* + R 2 ju . (11) 
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We assume the factorable dependence 5i = i?*f2*(i? — 
R*)C{t), where C(t) is a dimensionless function deter- 
mined subsequently. Consequently, 

8w(R,t) = n* R ~ R R * [C{t)-2] . (12) 

Therefore, equation (10) becomes 

or more compactly, 



Introducing 



d(k R AR) 
dt 



A R 

-SI.— [C(t)-2] 



Aw(t) = n,^ffi[r(t)-2] , 



allows us to rewrite this equation as 
d(k R AR) 



dt 



= -Alu 



(13) 



(14) 



(15) 



Because Ur oc — £?^(i?*) from equation (8), k R AR is pro- 
portional to the toroidal flux (in one direction) in the layer. 
Equation (15) expresses the fact that the rate of change of 
this flux is proportional to the 'rate of twisting,' Aoj(t). 

2.1. Angular Momentum Conservation 

The total angular momentum of the heated layer minus 
AMO,^ is simply 

f R,+AR m 

AL = RM* / dR{R - R*)£(t)— , 

Jr. dR 

= c M R*n*AMAR(t)£(t) . (16) 
Here, cm is a dimensionless constant of order unity assum- 
ing a self-similar mass distribution in the layer, dM/dR = 
f[(R-R*)/AR}. 

Changes in AL are due to twisting of the magnetic field. 
The vacuum field in the region R > R* + AR rotates but 
it is otherwise unaffected by processes within the heated 
layer. Thus the change in AL is due to the outflow of 
angular momentum from the surface of the star R = 
That is, 



dAL 

dt 



dS R Rsin0 



H(B«)R(B*) 



4ir 



R=R, 



^rRIBq , 



(17) 

where the fields are given by equation (8) and the integra- 
tion is over the area 2irR 2 of the equatorial region of the 
star. Using equation (16) gives 



dt^ = t^ 



(18) 



where 



AR 
~R 



1Z = — — , K, = kRR. t , 



R*Bq 



* Ac M AM ' 

where lor — const is an Alfvcn frequency of the heated 
layer. The actual oscillation frequency of a layer of con- 
stant thickness AR is larger than ujr by a factor of 
(R^/AR) 1 / 2 (see equation 21). 

For representative values from Cumming and Bildsten 
(2000), 



UjR 



3.16/s / B 



10 21 g 



1/2 



R* 



1/2 



(19) 



s/\~cm~ vio 8 g; " v AM J V 106cm / 

For a layer of constant thickness, the Alfven speed in the 
layer is va ~ (wb /tt) \/R* AR. For the given reference 
values and AR = 10 3 cm, va ~ 3.2 x 10 4 cm/s. 



2.2. Main Equations without Dissipation 

Using equation (14), equations (15) and (18) can be 
written as 



djUJC) 
dt 
dW 
dt 



£-211 



(20) 



where 



W = 



Aoj 
12^ 



is the dimensionless 'rate of twisting.' The initial condi- 
tions are that /C(0) = which corresponds to no initial 
toroidal magnetic field, and W(0) = —27^.(0) which cor- 
responds to the specific angular momentum of the layer 
being equal to its equilibrium value (i.e., £(0) = 0). 

Equations (20) are a linear system for (/C, W) in that 
lZ(t) is a given function. For 1Z — const, the solution of 
equations (20) is oscillatory with angular frequency 



^osc — 



LOB 

Vn 



(21) 



Thus the period of the oscillation is proportional to \/ AR. 
Equations (18) have the form of Hamilton's equations, 

Q= dH/dP, P= -dH/dQ, with P = w|/C£/ft* the 
canonical momentum, Q = W the canonical coordinate, 

and the Hamiltonian H = P 2 /{2TZ) + w|Q 2 / 2 - 2 ft P. 
These equations are idealized in the respect that they ne- 
glect the magnetic diffusivity and the viscosity of the layer. 
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Fig. 2. — The contours show dynamic power spectra of the X- 
ray burst of SAX J1808.4-3658 on October 18, 2000 (Figure 1 of 
Chakrabarty et al. 2003). The solid curves show the X-ray count 
rate. The bottom panel shows a close in view of the top panel. 



2.3. Equations with Magnetic Diffusivity 

For finite conductivity a of the layer plasma, the term 
?7 m V 2 B is added to the right-hand side of equation (4), 
where r\ m = c 2 /(4tt<j) is the magnetic diffusivity. In place 
of equation (6) we have 

= -(v'-V)* + ^A**, (22) 



dt 



where A* = d 2 /dR 2 + (1 / R 2 )d 2 / d<j) 2 is the adjoint Lapla- 
cian operator. 

We solve equation (22) with a generalization of equa- 
tion (7) where the radial wavenumber is complex, k R — > 
kft r + ikm. We introduce the dimensionless wavenumber 
{JC r ,JCi) = (k Rr R*, k Rl R*), and r m = R 2 /ii m - Thus 



d(WC r ) 

dt 
d{WCi) 



— — o*yv K, r K,i 



-(IC 2 + l-IC 2 ) , 



B 



2K 



(23) 



dt 
dW 
dt 

The initial conditions are that /C r (0) = 0, /Ci(0) = 0, and 
W(0) = -211(0). In contrast with equations (20), the 
evolution equations with dissipation are nonlinear. 



2.4. Magnetic Diffusivity 

Note that r m = R 2 /i] m is the Ohmic dissipation time- 
scale of the layer multiplied by (i?„ / AR) 2 . It is less than 
or equal to its value given by the classical Spitzer value of 



T m < TSpitz 



10 13 s 



i?* 



10 6 



cm 



/ ji \ 3 / 2 

{ Wk J 



(24) 



where T is the temperature of the heated layer, Z is an av- 
erage atomic charge, and the Coulomb logarithm has been 
taken to be 8 (Cumming & Bildsten 2000). 

Instabilities and associated turbulence in the heated 
layer are expected to give r m <C rs p itzor- Under differ- 
ent conditions the deep part of the layer may be con- 
vective and the top part radiative as discussed by Cum- 
ming and Bildsten (2000). The sound speed in the heated 
layer c s w 2.9 x 10 8 cm/s(T/10 9 K) 1 / 2 is generally much 
larger than the Alfven speed va estimated previously as 
va ~ 3.2 x 10 4 cm/s at the beginning of the burst. Later in 
the burst the magnetic field is twisted and the magnitude 
of the toroidal field increases substantially due mainly to 
the field being confined to the thin layer AR. Thus the 
Alfven crossing time tA = AR/va decreases from its initial 
value ~ 0.03s assuming AR — 10 3 cm. The strong toroidal 
field may lead to a buoyancy instability. Notice however 
that the layer has a strong shear with a velocity difference 
across it of order Av w 2£!„ AR ~ 5 x 10 6 cm/s so that the 
shearing time t s hcar = AR/ Av ~ 2x 10~ 4 s is shorter than 
tA in the initial part of a burst. Furthermore, notice that 
the upward buoyant motion of a blob will be strongly influ- 
enced by the Coriolis force 2pvx i2, where v is the velocity 
in the rotating reference frame. This force has a stabilizing 
affect on the blob motion and tends to give circular 'gyro' 
motion of radius r g = |v|/2£l ss 6 cm [|v|/(3 x 10 4 cm/s)] 
about the fl direction (Tritton 1988) for Q,*/2tt = 400 Hz. 
Thus, r g for Alfven speed motions is smaller than the layer 
thickness AR even for conditions where the field strongly 
increased. 

Because da C c s , the layer may also be unstable to the 
magnetorotational instability (Balbus & Hawley 1998; Ve- 
likhov 1959; Chandrasekhar 1960; Moiseenko, Bisnovatyi- 
Kogan, & Ardeljan 2006). The saturation of the insta- 
bilities is assumed to give rise to a turbulent diffusivity 
T] m t much larger than the Spitzer value and a turbulent 
viscosity v t ~ 7] m t (Bisnovatyi-Kogan & Ruzmaikin 1974). 
The viscosity is taken into account in the next subsection. 
The turbulent viscosity and diffusivity can be written in 
a form analogous to the Shakura-Sunyacv (1973) "alpha" 
prescription for a thin disk taking into account that the 
maximum eddy size is AR and the maximum buoyant ve- 
locity is ~ va with va the initial Alfven speed. Later in 
a burst, va increases substantially but at the same time 
the mentioned effects of shear and Coriolis force act to 
limit the motion of buoyant blobs. Therefore, we estimate 
Vmt — cuvaAR with a < 1 which gives 



3 x 10 7 - 



VA 



AR 



10 3 cm / V 3 x 10 4 cm/s 



(25) 

Note that for r) m = 3 x 10 7 cm 2 /s, the magnetic diffusion 
time across the layer is tdiff = AR 2 /r/ m ~ 0.033s is of the 
order of the initial Alfven crossing time. 
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We emphasize that the values of r) mt and v t are highly 
uncertain in that relevant linear and nonlinear theory and 
simulations have yet to be done. An extensive literature 
exists on the the theory and simulations of the related 
problem of the stability and motion of magnetic flux tubes 
in the convection zone and overshoot region of the Sun 
(e.g., Schusssler et al. 1994). It appears possible that 
these methods can be adapted to the magnetic field sta- 
bility of the X-ray burst sources. 

2.5. Equations with Diffusivity and Viscosity 

Viscosity of the heated layer is taken into account by in- 
cluding on the right-hand side of equation (17) the viscous 
angular momentum flux across the surface R = R*. This 
term is 

dS R R sin 9T^, = - 2-nRipv ^ 



L dR lR=R, 

-2nRip(R*) 



(26) 



where Tj^ is the viscous contribution to the momentum 
flux-density tensor. Thus, the equations for the layer in- 
cluding Ohmic and viscous dissipation are 



d(WC r ) 

Jt 
d{HKj) 
dt 
dW 
dt 



-n,yv K-rK-i , 

±(JC 2 r + l-JC 2 ) 

Tm 



1_ vv 

T„1¥ 



2K 



where 



where 



Rl 



(27) 



(28) 



= cmAM/[2ttR 2 ARp(R^)] is a dimensionless 
constant of order unity if the mass distribution of the disk 
is self-similar. Numerical solutions of equations (26) indi- 
cate that the influence of viscosity is negligible for r v >r m . 

3. SAMPLE SOLUTIONS 

To solve equations (27) we need 1Z(t) = AR(t) / i?* which 
is not known from observations. A rough estimation of this 
function can be made assuming AR oc T and the X-ray 
intensity I oc T 4 . The X-ray intensity of the SAX J1808.4 
burst shown in Figure 2 (Chakrabarty et al. 2003) is fit- 
ted approximately by / oc 1/{1 + [(t — 5)/6.5] 2 } and for 
10 < t < 40 accurately by I oc 1/t 2 . The corresponding 
layer thickness is TL = ft(0)/{l + [(t - 5)/6.5] 2 } 1 / 4 , with t 
in seconds. 

Figure 3 shows a sample case relevant to the observed 
burst shown in Figure 2. Figure 4 shows the trailing spi- 
ral magnetic field at t = 0.4s for this case. Later, when 
fC r < 0, a leading spiral field forms. The ratio (3 = B^/Br 
at R = i?* is (3 = —K, r and this is seen to vary from —1450 
at 0.42s to 420 at 1.45s. The first positive peak of Aw/fi» 
occurs at t\ w 0.25/ujb for small damping, r m > 10 5 s. 
The time between the first and second positive peaks of 
Auj/VL* is Ati2 w 0.53/w B . 

Note that r m = 10 5 s corresponds to a diffusivity r\ mt — 
Rl/ T m = 10 7 cm 2 /s for i?* = 10 6 cm. From Figure 3 it 
is seen that the twisting amplifies the toroidal magnetic 
field to a peak value of the order JC r sa 1400 larger than 



the initial poloidal field of the star. The amplification is 
due to the field being confined to the thin layer AR rather 
than it being wrapped many times around the star. For 
the case of Figure 3, the field is wrapped two turns in the 
clockwise direction during the first ss 0.7s, and then it is 
unwrapped, wrapped, etc., subsequently. The magnetic 
diffusion time across the layer is idiff = AR 2 /r] m w 5.6s. 

Figure 4 shows the "phase slip," J * dt'u>(t') — £l*t, dur- 
ing a burst for the same conditions as Figure 3. Figure 5 
shows the nature of the wrapped field at t = 0.4 s for the 
case of Figure 3. 

Figure 6 shows Aw/fi* over a longer time interval for 
sample cases. The period of the oscillations of Auj is pro- 
portional to \/AR which in this case is proportional to 
1/t 1 / 4 . 



0.02 



0.01 




-0.01 



-0.02 



Fig. 3. — Sample solution of equations (23) for W(t) = Aw/H*, 
Kr(t) and /Q(t) relevant to the burst shown in Figure 2. The con- 
ditions arc u)b = 0.252/s and r m = 10 5 s, which corresponds to a 
magnetic diffusion time across the layer of t^g = Afi 2 /ij ra t ~ 5.6s. 
Also, TZ(t) is given in the text with TZ(0) = 0.0075 and Au(0)/fi, = 
0.015. 
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Fig. 4. — The "phase-slip during a burst for the conditions of 
Figure 3, where <f>(t) = f* dt'w{t'). 
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Fig. 6. — Sample solution of equations (23) for W(i) = Au/fl, 
for a longer interval for (a) u)b = 0.0504/s and r m = 10 6 s, (b) 
uj b = 0.252/s, r m = 10 6 s, and (c) w B = 0.252/s, r m = 10 5 s, with 
other conditions the same as for Figure 3. Note that T m = 10 5 s 
corresponds to a magnetic diffusion time across the layer of t^ifi = 
AR 2 /r] m t ~ 5.6s. For r m = 10 6 s the diffusion time is ten times 
longer. The period of the oscillations of Au> decreases gradually 
with time. 



4. CONCLUSIONS 



50AR 




Fig. 5. — Field lines in the heated layer at t = OAs for same 
case as Figure 3. The radial thickness of the layer has been 
expanded by a factor of 50. The field line is given by R*<f>' = 
-K r (R - R*) - (Ki/K r )R* ln{abs[cos(/C r (i? - R*)/R*)}}, where 
IC r = 1447, K t = 597., and AR/R t = 0.00678. N and S indicate 
the north and south magnetic poles. 



We have derived a simple model for the influence of a 
neutron star's magnetic field on the rotation of its surface 
layer rapidly heated by thermonuclear burning. The model 
assumes that the star's magnetic moment is perpendicular 
to its rotation axis. The burning causes the expansion of 
a thin outer layer of the star, AR(t). The layer rotates 
slower than the star due to angular momentum conserva- 
tion. The shear between the star and the layer acts to 
twist the star's magnetic field giving at first a trailing spi- 
ral field. The twist of the field acts in turn to 'torque up' 
the layer increasing its specific angular momentum. As 
the layer cools and contracts, its excess specific angular 
momentum causes it to rotate faster than the star which 
gives a leading spiral magnetic field. The oscillation of the 
angular velocity of the layer is a result of the tension of the 
twisted magnetic field. Non-uniformity of the star's photo- 
sphere (at the top of the heated layer) is due the magnetic 
field and this gives rise to the observed X-ray oscillations. 
The fact that the layer periodically rotates faster than 
the star means that the X-ray oscillation frequency may 
"overshoot" the star's rotation frequency. Observations 
by Chakrabarty et al. (2003) of an X-ray burst of SAX 
J1804. 4-3658 show clear evidence of the overshoot of the 
frequency of the X-ray oscillations. 

Equations of the model are for the difference in angular 
velocity between the layer and the star and the radial and 
toroidal components of the magnetic field. The frequency 
of the oscillations is proportional to the initial poloidal 
magnetic field of the star, inversely proportional to the 
square root of the mass of the heated layer, and inversely 
proportional to the square root of the layer's thickness. In 
the absence of magnetic diffusivity and viscosity, the equa- 
tions are linear and constitute a Hamiltonian system. The 
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equations become nonlinear with the magnetic diffusivity 
included, but the inclusion of viscosity adds a linear term. 
The diffusivity and viscosity are probably due to turbu- 
lence in the heated layer but the level of the turbulence 
is highly uncertain. The model has two important param- 
eters, one is the oscillation frequency proportional to the 
initial magnetic field, and the other is the damping time 
due to magnetic diffusivity. The value of the magnetic dif- 
fusivity is highly uncertain. We find that the twisting can 
amplify the toroidal magnetic field to a peak value of the 



order 10 3 larger than the initial poloidal field of the star. 
The amplification is due to the field being confined to the 
thin layer AR rather than it being wrapped many times 
around the star. 
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